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Abstract 

Reported in this paper is the impact of the fluctuations of the ge- 
ometry of the nano-meter gas containers in the medium on the NMR 
line shape of the gas inside of the nano-containers. We calculate ex- 
actly the NMR line shape of the gas of spin-| carrying molecules for 
two typical dynamics of the nano-container volume and the orienta- 
tion with respect to the external magnetic field: (i) for a Gaussian 
stochastic dynamics, and (ii) for the regular harmonic vibrations. For 
the Gaussian ensemble of static disordered containers having an infi- 
nite correlation time, t c — > oo, the overall line shape is shown to obey 
a logarithmic low frequency asymptotics, I(oj) = const x hi(jj), at 
u — > 0, and exponentially decaying asymptotics in a high frequency 
domain. For the Gaussian ensemble of the rapidly fluctuating con- 
tainers of a finite r c , the overall line shape has a bell-shaped profile 
with ~ lu~ 4 far wing behaviour. In addition, we calculate exactly a 
satellite structure of the NMR line shape when the nano-bubbles in a 
liquid are affected by the harmonic deformations due to the acoustic 
waves. 
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1 Introduction 



The systems we are studying in this paper are of interest due to the recent 
experiments [T in which NMR responses from the H2 gas confined in the 
nano-meter porous hydrogenated silicon have been used for measuring the 
size and the orientation of the nano-scale pores The accessibility of 
these measurements is relied on the explicit dependence of the line width of 
the NMR line shape on the confining volume and on the orientation of the 
ordered pores with respect to the strong external magnetic field. When a 
gas of spin-carrying molecules is confined within a nano-meter sized region, 
the motionally averaged (effective) dipolar interactions between the nuclear 
spins take on the finite values due to the spatial averaging of the dipolar 
interactions over the finite ( nano-meter sized ) region, hence, allowing the 
estimation of the line width pQ by the Van-Fleck formula. In addition, the 
motionally averaged (effective) dipolar interactions lead immediately to the 
exactly solvable effective spin Hamiltonian (see Sec. 121 below) allowing the 
exact calculation of time dependent NMR responses. In our previous paper 
[3] we calculated the time course of the longitudinal polarization of a gas 
of nuclear spin-| carrying molecules confined in the static nano- containers. 
The purpose of the current paper is to extend the developed formalism for 
predictions of the NMR line shape in the presence of the dynamical nano- 
containers. The example systems we are exploring here are the cavitation 
bubbles that could be produced in experiments like those in sonoluminescence 
[3], the gas vesicles emerging in the non-invasive measurements of the liquid 
flow by the NMR tomography [3], 0, the bubbles in blood [7], the gas 
bubbles under the nucleation in the course of liquid-to-gas phase transitions 

nano-sized free- volume holes of polymers, rubber, etc. under the elastic 
deformations, the gas containers within the vibrating nano-tubes 0, |10j . 

It is common that the NMR signals offer the spectroscopic way of the 
measuring of the characteristic sizes and the relaxation times of the container 
that surrounds and traps the nuclear spins jllj . JT2*] . In studying the NMR 
signals from the iV spin ensemble with no special symmetry, the exponential 
with the number N of allowed unknown basis set of the spin Hamiltonian is 
to be incorporated to describe the NMR signals completely. However, the 
complexity of the spin Hamiltonian is greatly reduced when the nuclei are 
subjected to the fast thermal motion within the nano-meter container. Under 
these conditions it arises a vast difference of the NMR time scale t nmr = 10~ 4 
- 10 -5 sec, that characterizes the flip-flop transitions at the nearest proton- 
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proton distances, and of the motional time scale, t m = 1CT 11 - 10~ 12 sec, 
relevant for the round trip of the hydrogen molecules within the nano-meter 
container at room temperature. The presence of a reliable small parameter 

e = t m /t nmr « 1(T 7 (1) 

allows to average the underlying dipolar Hamiltonian over the spin spatial 
coordinates on the coarse-grained time intervals At obeying the inequality 

t m < At < i nmr , (2) 

giving rise to the motionally averaged spin Hamiltonian H ( see Eq-s (J3J and 
(jU) of Sec. |2J) with a unique spacing independent effective dipolar coupling 
PP, The exact spectrum of the motionally averaged spin Hamiltonian H 
has led to the development of the methods for characterizing the exact NMR 
spin dynamics, especially the non-ergodic spin dynamics and the line shape 

We mention also the recent studies ^3], 03j in the area of the NMR re- 
sponses from a confined gas in the framework of the phenomenological Bloch- 
Torrey equation as well as the NMR measurements |T5|, ^H], \U\ revealing 
the reduction of the atomic mobility with respect to the bulk mobility due 
to the confinement effect. Although the previous papers have been largely 
concerned with the NMR measurements in the static nano-containers, the 
methods developed can be taken over for the cases of "flexible" walls of 
nano-containers. The point is that the dynamics of the surface of a typical 
bubbling behavior in liquids occurs at a millisecond time scale i.e. at the 
same time scale which is relevant for the NMR spectroscopy. In Sectional we 
give a general formalism of the free induction decay (FID) for an arbitrary 
moving nano-containers. This is followed by analysis in Section El of the FID 
and the line shape for the Gaussian temporary fluctuations of the volume 
and of the orientation of the nano-containers. Our concern in the Section |H] 
is the line shape for a wide range of the correlation times and the amplitudes 
of the fluctuations of the nano-containers. Section Ogives the line shape from 
the nano-containers subjected to a regular harmonic vibrations at a single 
frequency as well as the line shape from the nano-containers subjected to 
the harmonic vibrations with the Gaussian distributions of the frequencies. 
Finally, Section El summarizes the major conclusions of the calculations. 
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2 Line shape of a gas within nano- containers 
with a time-dependent volume 



Consider a gas of N spin-^ carrying molecules confined in a moving nano- 
containers in the strong external magnetic field B. On the coarse-grained 
time intervals At (0) the effective spin dynamics is described by the motion- 
ally averaged spin Hamiltonian , 

H = u I z + ^D(t)(3I 2 z -I 2 ), (3) 

where the u = 7 p I? with 7 p standing for the proton gyromagnetic ratio, 
the nuclear spins are specified by the spin-| operators I ia , i = 1,...,N, 
(a = x,y,z), the operators I a = Y^f^ha-, are referred to as the projections 
of the total spin operator onto the x, y, z axes, respectively. Next, we assigned 
I 2 = I 2 + I 2 + I 2 to the square of the total nuclear spin operator. Finally, 

the motionally average (effective) dipolar coupling between all ( pairs of 
the spins in the nano- container is 

D(t) = ft^(3<x*?e(t)-l). (4) 

Here, the V(t) is the volume of the nano- container, the dimensionless form- 
factor f(t) depends on the shape of the nano-container and 9(t) denotes the 
time dependent orientation of the nano-container with respect to the external 
magnetic field B, see Fig. 1. The equivalence of the effective coupling D 
for all the pairs of the nuclear spins is due to the fact that all the nuclei 
inside the nano-container are involved in the equivalent fast thermal motion 
inside the nano-container over the long NMR time scale t nmr . In the absence 
of the nano-container's motion, the nuclear motion inside the static nano- 
container gives rise to the unique motionally averaged spin coupling D pQ, 
jsJj. The analytical dependence of the coupling D on the nano-container's 
volume V admits an immediate extension of the coupling D to the case of 
the time dependent volume V by invoking the adiabatic framework 0, i.e. 
by regarding the motion of the nano-container to be slow as compared to the 
fast thermal motion of the nuclei within the nano-container. 

The explicit form of the function f(t) for the ellipsoidal container is in 
Ref. [Sj. For the static containers of the nano-meter volume, V ~ 10 3 A 3 , 
the effective coupling in Eq. (@J) is evaluated as 10 3 times smaller than the 
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characteristic flip-flop coupling (~ t n ^ r = j^h/a^ ~ 10 4 — 10 5 Hz ) of two 
protons at a nearest separation ao ~ 1A. This reduction of the dipolar 
coupling is referred to as the fast motionally narrowing (in 10 3 times) of the 
line shape as compared to the line shape of the static nearest protons. 
The line shape is the Fourier transform of the FID [llj, 



F(t)=tr{p(t)I-}/tr {/+/_}, (5) 

where I± = I x ± il y and p(0) = 1+ is the initial density matrix in the 
high temperature approximation |18j . The density matrix p(t) of the iV-spin 
ensemble in the rotating reference frame obeys the Liouville-von Neumann 
(L-vN) equation (% — 1) 

(^[^(aij-np] (6) 

In solving L-vN equation ©, we introduce the phase shift, 

<p{t) = \f dt'D(f), (7) 



2 Jo 

then, the L-vN equation © is solved to be 

p{t) = e MW I+e -MW = e <3 V (t)(2J,-i) /+- (g) 



In deriving (jHJ), we put 7i = 3I| — I 2 , use the commutators 
the identity 



/ I 2 



and 



jWl I+e -e*pi> = e ^(2/,-i) /+ _ (g) 

The trace in Eq. (jSJ) is easily performed in the total occupancy number 
representation of N\/(N^ !iVjJ)-fold degenerate basis set |A^-, A^), where A" = 
A^ + Aj and 7 Z = |(A"-f — AT|). The trace gives the sought FID for an arbitrary 
time-dependent coupling D(t), 

F(t) = (cos^t)))^ 1 . (10) 



On the NMR reasonable time scale t < t nmr ~ 10 4 sec, we get (D)t < 10 2 , 
so f(t) <C 1 and the FID, F(t), of Eq. (fTU|) transforms into 



PQ) _ e (iV-l)ln(co S (3 V (t))) ^ e -f (3v>(t))- /-q) 
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for a large number of spins, N 3> 1, in the nano-container. The FID, F(t), 
of Eq. (jlljl involves the effective coupling D(t) of Eq. (@J) as the input 
parameter to the phase shift ip(t) of Eq. (J2J. Varying of the function D(t) 
yields a variety of the models of the NMR line shape which are commonly 
discriminated into the two major groups, viz. the models of the homogeneous 
or inhomogeneous line width ^T], ^2]- in the following, we explore the FID 
for two dynamical scenarios of the container motion, for a stochastic Gaussian 
dynamics (Section |3J) and for the regular harmonic oscillations (Section H] ) . 



3 Line shape from fluctuating nano-containers 

When the nano-sized containers are sensitive to the fluctuations of the envi- 
ronment, we are free to calculate the line shape by assuming the Gaussian 
fluctuations of the coupling D(t) (jU, 

D(t) = (D) + 5D(t), (12) 

with the SD(t) standing for the Gaussian random noise characterized by the 
first two moments 

(8D{t)) = 0, (SD^SDfa)) = ((5D) 2 ) C{\t x - t 2 \), (13) 

where ((5D) 2 ) is the variance of the fluctuations and the C(t) denotes the 
correlation function, for example, C(t) = exp(— t/r c ), with r c being the cor- 
relation time. The averaging of the function F(t) over the Gaussian fluctu- 
ations SD(t) is carried in two steps: first, we rewrite the FID, F(t),ofEq. 
(fTTj) by introducing the Gaussian parameterization, 

= r+°° dx ^.wsvvfl (M) 

J-oo \/7l 



this is followed by the second step of applying the formula for the averaging, 
see e.g. jTI], of the function F(t) of Eq. (JUjl over the random Gaussian 
process SD(t) entering the phase (p(t) of Eq. (j7J), 



exp (-in J* 5D(t')dt'^ = exp(-K 2 (((5 J D) 2 )r 2 (t)^ 

T 2 (t) = (\t-t')C{t')dt\ (15) 
Jo 
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with the constant k = 3xyiV/2. Averaging by Eq-s. (|T5j) . (fTlj) gives the 
sought FID 



As it stands, the FID, F(t), of Eq. (jlfij) encodes an information on the mean 
volume and the mean orientation of the cavity with respect to the external 
magnetic field as well as on the fluctuations of the cavity. 

If the fluctuations were absent (a = in Eq. (|16j) ) then the FID, F(t), 
of Eq. (jTfi)l forms the line shape, 

1 r 00 

I(uj) = - F(t)cosutdt (17) 



7T JO 



of the Gaussian type, I(uS) = j^expy — j^J. Fluctuations of the cavity 
(a 7^ in Eq. (jTfij)). result in the broadening of the line shape so that the 
second moment, see e.g. [TT], reads 



d(tf 



9/V r i 

— \{Df + {{5Df) . (18) 



t=o 



In deriving the second moment, use is made of the expression T 2 in Eq. 
and the property C(0) = 1. 

To make the calculations of the line shape I {(J) more explicit, we take 
the exponential correlation function, C(t) = exp(— t/r c ). It follows 

T 2 (t) = r c 2 (expH/r c ) + --l). (19) 

In the analysis of the FID, F(t), of Eq. HBJ) with the function T 2 (t) of Eq. 
(HHJ) , it seems very useful to consider the temporary fluctuations on two very 
different time scales. 

When r 2 ((5D) 2 ) < 1, the function T 2 (t) is T 2 {t) = r c t for times t > r c , 
so that the FID of Eq. dEJ) reads F(t) = exp(-^-t) /VI + ar c z/ 2 t, and its 
Fourier transform of Eq. (JT7j) admits the following representation 

1 € z 

I(u) = — = — - Ref -= erfc(Vi) J , z = (2ar c u)- 2 + za;(ar c ^ 2 ) -1 , (20) 
aV 7rr c l/ v 2 
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in terms of the function erfc The function I(oS) of Eq. (|2U|) has a bell- 
shaped profile with an intermediate Lorentzian asymptotics I{ui) = -r/(T 2 + 
uj 2 ) where T = (4ar c ) _1 ( see in [TH] the asymptotics of the function erf c 
at \y/z\ > 1 or r c ^ ). Far wing calculations of the line shape I(u) require 
the FID F(t) at <— t < r c , which is provided by the function T 2 (t) of 
Eq. (jTUj) at t — > 0. The sought asymptotics of the line shape I(u) can be 
evaluated as 



Q/V 1 



UJ 



oo 



(21) 



on integrating the line shape I{uj) of Eq. (JT7j) four times by parts and em- 
ploying the derivatives F'(0) = 0, C"(0) = — r" 1 . 

On the other hand, when t 2 ((5D) 2 ) ^> 1, the function T 2 (t) becomes 
T 2 (t) = t 2 /2 for < t < t c -y oo. The function T 2 (t) leads imme- 
diately to the slowing down of the FID of Eq. (j!6|) in the form F(t) = 
exp(— ^\ / \Jl + \avH 2 . Consequently, one is left with the line shape 




00 cos ut dt 



\au 2 t 2 



ITU 



2 T ^ f ujV2 



and 



e 2 Q 1 / uj\pl 
expl 



CO 



7iv\/2a 



UJ 



V\ OL 



UJ — ► oo. 



(22) 



(23) 



In deriving the I(uj) of Eq.-s (|22j) and ()23|1. use is made of the integral repre- 
sentation of the modified Bessel function K (x) and its asymptotics ^HI- in 
order to shed light on the slowing down of the FID F(t) at t — > oo (and, thus, 
on the logarithmic singularity of the line shape I(uj) at uj — > in eq. (|22jl ). 
we regard a disordered almost static (frozen at r c — > oo) distribution of the 
nano-containers with the various volumes. The FID from an individual nano- 
container is described by the Gaussian-in-time function with the relaxation 
rate v = 3D^jN/2. The main contribution to the sum of the individual FID's 

weighted by the Gaussian probability density (2Tr(D 2 ))~ 1 ^ 2 e~ D2 / 2l ^ D2S} of the 
static fluctuations of the coupling D comes from the fluctuations having the 
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D = 0, i.e. from the containers with a large volume or oriented at the magic 
angle arccos(l/v / 3) ( see the expression for the coupling D of Eq. (J3J)). The 
FIDs from the large cavities are slowly damping providing a slow damping 
of the overall signal F(t) ~l/iatt^oo rather than the Gaussian-in-time 
asymptotics. Notice that the slowing down of the overall FID from the Gaus- 
sian ensemble of the static fluctuating nano-containers akin the slowing down 
of the unimolecular decay on the static disordered traps [201, |2Tj - Fig. 2 
is aimed to show the appearance of a low frequency singularity of the line 
shape while increasing r c to infinity. The line shape for fixed t c u = 10 at the 
various a is displayed on Fig. 3. In addition to having the broad shape, the 
line shape shows the singularity at the zero frequency asymptotics at 
the large fluctuations a. 

It is worth to gather a small number of characteristic quantities involved 
in the Figure 3; for N = 500 spin-| molecules ( j^h = 2tt ■ 120 Hz • nm 3 ) 
within V = 45 nm 3 nano-container and the form-factor f(t) ~ 2 of Eq. (£Q), 
the motionally average dipolar interaction is evaluated as (D) = 2-k ■ 5.3 Hz, 
thus, v = 0.25 • 10 3 Hz and r c = z/" 1 ~ 4 • 10" 3 s. 

The lessons drawn from this Section are that the fluctuations of the nano- 
containers give rise to the deviation of the line shape from the standard 
Gaussian and the Lorentzian shapes, and that for almost static disordered 
nano-containers at t 2 ((5D) 2 ) ^> 1, the line shape I{oj) gets narrower at uj — > 
and broader at uj —>■ oo as compared to the bell-shaped profile occurring 
for frequently fluctuating nano-containers at t^((5D) 2 ) 1 . 

4 Line shape from vibrating nano-containers 

Acoustic waves in the liquid surrounding the nano-bubbles can induce a syn- 
chronized harmonic vibrations of the nano-bubble volumes and their orien- 
tations 0], thus, affect the NMR line shape if the acoustic waves are at the 
NMR relevant frequency domain 1 — 10 kHz. The same physical picture of 
the NMR responses should appear for a gas within the vibrating nano-tubes 
0, [ID]- For all these vibrating nano-containers, we can regard the coupling 
D(t) of Eq. (@J) to be a harmonic function of the time, 

D(t) = (D)(l + ecos(nt)) (24) 

with parameter e < 1 assuming a weak vibrations of the nano-bubble volumes 
and the orientations. On the NMR time scale t ~ 10 -3 sec, the phase shift 
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(p(t) of Eq. (J7J) becomes 

v(t) = l{D)(t + ^sm(nt)), (25) 

so that yj(t) <C 1 and the signal F(t) of Eq. (jTQ^ again transforms into the 
F(t) = e~2( 3Ar ^(*)) 5 ( see Eq. (IlljO . For a weak vibrations, it is enough to 
expand the function F(t) in the powers of the parameter e keeping only the 
terms up to e 2 and linear in a small factor v 2 = 1(D) 2 TV, ( by Eq. (gj), the 
coupling (D) ~ 1/iV), giving 

^)=e-^ /4 (l-^sin(^)-, 2 Q 2 sin 2 (^)). (26) 

The Fourier transformation brings the signal F(t) of Eq. (|2l)j) into the line 
shape 

i(uj,n) = -^ e ~^+ ]T (/i(w,n) + /i(-w,n) 



= -wf( 1 + o) 



^ Wife) (27) 

In the absence of the bubble vibrations (e = 0), the line shape I(u>, Q) exhibits 
the motionally narrowed peak at frequency uj = with the line width 2u. 
Weak bubble vibrations at the single frequency Q leads to the appearance of 
the symmetric satellite pairs of the line shape I(u,Q) (|27|) at the multiple 
frequencies uj = ±Q, ±2Q. By accounting for the n-th term in the powers 
of the amplitude e, the satellite pairs at frequencies uj = ±nQ, n — 3, 4, . . . 
arise. 

Now, let an ensemble of many individual gas bubbles is spread over the 
liquid having a random local vibrational frequency Q due to an intimate 
fluctuations of the liquid. Under these conditions, the nuclei belonging to 
different bubbles are (indirectly) subjected to a different local vibrational 
frequency resulting to the inhomogeneous broadening of the NMR spectrum. 
We think of the ensemble of the bubbles as a single bubble affected by the 
vibrations with a continuous Gaussian distribution of frequencies 

(n-n ) 2 I 

D(Q) = A Q 2 e sr- with A = -. r (28) 
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ensuring the normalization /f^ dQD(Q) = 1. The pre-factor f2 2 in Eq. (|28|) 
is taken for ease of performing the averaging of the line shape I(u, Q) of Eq. 
(|2"Tj) over the distribution D(Q) of Eq. 



yielding 



/oo 
dQD(Q)I(u,fl,), (29) 
-oo 



(/(«)> = ^=e"^+ + 



i=l,2 



2 1/(1 + 5 2 ) 3 / 2 L2 ' A 2 (l + <5 2 )""" 



("+ n o> 2 

e ^Ta 2 ; 



e 2 Aiz/A t^+2n ) 2 
21 j 16(l + 5 2 ) 1 /2 e ' ^ 

with 5 = A/V. The overall line shape (/(&>)) of Eq. is still narrow at uj = 
and reveals the two symmetric satellite pairs at frequencies ui = ±f2o, ±2fio 
having a broad line width 2a/ z/ 2 + A 2 and 2V i/ 2 + 4A 2 , respectively, see Fig. 
4. 

In general, the satellite pair of the amplitude e n is described (to within the 



pre-exponential factor) by the Gaussian shape with the line width 2a/ v 2 + n 2 A 2 
2nA, i.e. 2n times larger than the dispersion, A, of the frequencies in the 
spectral density D(Q) of Eq. (J2%Jl . Thus, the position of the n-th satellite 
pair at u = ±nf2o and the broadening of the n-th satellite, 2nA, provide 
the NMR spectroscopic characterization of the nano-bubble vibrations hap- 
pening at the mean vibrational frequency Qq and with the dispersion of the 
frequencies A. 



5 Conclusion 

The focus in the paper is on the exact NMR line shape theory of a gas of 
spin-1/2 carrying molecules confined within the fluctuating nano- containers. 
Two typical dynamics of the nano-containers was treated, viz. the Gaussian 
stochastic dynamics and the regular harmonic vibrations. 

(1). Of the variety of the Gaussian random fluctuations of the nano- 
containers, the most striking fluctuation effect on the NMR line shape is due 
to the fluctuations at the large correlation times, r c — > oo, (for almost frozen 
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disordered ensemble of the various nano-containers) and at the large ampli- 
tudes of the fluctuations of the volume and orientation of the nano-containers. 
Under these conditions, the NMR line shape behaves as I{uj) = const x In - 
at uj — ► and exponentially decaying at the large frequencies, uj — > oo. Al- 
ternatively, when the conditions are specified by the small correlation times, 
t c — ► 0, or at the small amplitudes, a, of the Gaussian fluctuations of the 
nano-containers, then the line shape has the bell-shaped profile with the 
power law ~ uj~ 4 at far wings. The line width and its precise shape specify 
the mean volume, the mean orientation of the cavities as well as the deviation 
of the volumes and the orientations from the mean values. 

(2). If the driving sources of the vibrations support the harmonic vibra- 
tions of the bubble volumes and of the orientations at a single frequency Q, 
then the line shape has the spike satellite pairs with a narrow line width 2v 
at the frequencies uj = ±nO, n — 1, 2, . . . around the central spike at uj = 0. 
For the Gaussian distribution of driving frequencies with the mean Qq and 
the dispersion A, the central spike at uj = remains to be narrowed with 
the line width 2u, however, the satellite pairs at uj — ±nf2 , are subjected to 
broadening in the way that the n-th pair has the line width ~ 2nA. 

The upshot is that the paper demonstrates how the fluctuation dynamics 
in the medium can be characterized by the NMR spectroscopy of the gas 
within the fluctuating nano-containers. 
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Captions to figures. 

Fig. 1. Cartoon of the two positions of the nano-container which is mov- 
ing in a liquid undergoing the deformations of the volume and/or the varia- 
tion of the orientation 9{t) with the time. The nano-container confines the 
gas of nuclear spin-| molecules uniformly spread inside the nano-container. 

Fig. 2. The line shapes of Eq. (fTTj) are calculated by the Fourier 
transform of the FID F(t) of Eq. (JTED with the function T 2 (t) of Eq. (JTHD for 
the various values of the combination t c u. The parameters v and the a = 1 
are from Eq. (fTBj) . 

Fig. 3. The line shapes I(u) of Eq. £□$ for the FID F(t) of Eq. (fToT) 
with the function T 2 (t) of Eq. at the fixed r c z/ = 10, but the amplitude 
of the fluctuations are allowed to vary from a = 0.01 to a = 100. 

Fig. 4. The absolute value of the homogeneous (dashed-dot) and inho- 
mogeneous (solid) NMR line shapes of the vibrating bubbles with the nuclear 
spin- 1 molecules inside. The absolute value of the inhomogeneous line shape 
\ from Eq. (J3TJ|) is shown for the parameters e = 5, A = 2, Q = 2n, 
with all the frequencies being in the units of the frequency v of Eq. (|16|). 
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